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Ultra-peripheral collisions (UPC) of heavy-ions involve long range electromagnetic in- 
teractions at impact parameters twice larger than the nuclear radius, where no nucleon- 
nucleon collisions occur. The first measurement of photoproduction of J/ip and of 
two-photon production of high-mass e + e~ pairs in ultra-peripheral nucleus-nucleus in- 
teractions will be presented, using Au+Au data at yJs NN — 200 GeV. The measured 
cross sections at midrapidity are consistent with various theoretical predictions. [I] 

1 Introduction 

The study of photoproduction at hadron colliders has attracted an increased interest in 
recent years [2j EJ H] • Electromagnetic interactions can be studied without background from 
hadronic processes in ultra-peripheral collisions without nuclear overlap. This study focuses 
on the measurement of exclusively produced high-mass e + e _ -pairs in Au+Au collisions at 
y/s NN = 200 GeV, Au + Au — > Au + Au + e + e~ at midrapidity. The results [5] have been 
obtained with the PHENIX detector [5] at the BNL Relativistic Heavy Ion Collider (RHIC). 

The electromagnetic field of a relativistic particle can be represented by a spectrum of 
equivalent photons. This is the Weizsackcr- Williams approximation. The virtualitics of the 
equivalent photons when the field couples coherently to the entire nucleus are restricted by 
the nuclear form factor to Q 2 = (cu 2 /j 2 + q 2 A < h/R\, where 7 is the Lorentz factor of the 
beam and Ra is the nuclear radius. At RHIC energies, 7 = 108 and the maximum photon 
energy in the center-of-mass (lab) system is u) max ~ 3 GeV corresponding to a maximum 
photon-nucleon center-of-mass energy of W™ff x « 34 GeV. 

The exclusive production of an e + e _ -pair can proceed either through a purely electro- 
magnetic process or through coherent photonuclear production of a vector meson, which 
decays into a dilepton pair. The Feynman diagrams for the two leading order processes are 
shown in Fig. [TJ 

The production cross section of vector mesons is a good probe of the nuclear gluon 
distribution, Ga(x, Q 2 ), as well as of vector-meson dynamics in nuclear matter [7HHJI2]. For 
J/V'-production, PHENIX acceptance at y = corresponds to a mean photon energy of 
(Ej) = 300 GeV and nuclear Bjorken-x values of x — , m 2 j^jW 2 A ~ 1.5 • 10~ 2 . 

Measurements of coherent photonuclear production of the p meson [10], as well as 77 
production of low-mass pairs |11| have been performed by the STAR collaboration. The 
PHENIX analysis presented in Rcfs. [T^l H] and summarized in this paper is the first on 
heavy final states in ultra-peripheral nucleus-nucleus collisions. The cross section for J/ip 
and e + e~ photoproduction are compared with model calculations [Ul IT31 RBI HfJl 1TB] . 

2 Experimental analysis 

The data presented here were collected with the PHENIX detector at RHIC during the 
2004 Au+Au run at ^/s^7 = 200 GeV (Run-4). The PHENIX detector is equipped with 
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Figure 1: Lowest order Feynman diagrams for exclusive photoproduction of J/i[> (left) and 
dielectrons (right) in ultra-peripheral Au+Au collisions. The photons to the right of the 
dashed line are soft photons that may excite the nuclei but do not lead to particle production. 



multi-layer drift chambers (DC) followed by multi-wire proportional chambers (PC) with 
pixel-pad readout. The tracking arms also contain Ring-Imaging-Cerenkov (RICH) detectors 
and electromagnetic calorimeters (EMCal) for electron and positron identification. 

The events used in this analysis were collected by a level- 1 Ultra-Peripheral Collision 
(UPC) trigger set up for the first time in PHENIX in Run-4 as follows. A veto on coincident 
signals in both Beam-Beam Counters (BBC) selected exclusive-type events characterised by 
a large rapidity gap on cither side of the central arm. The EMCal and the RICH decctor 
were used to form a trigger (ERT) to select events with at least one of the two high-energy 
e ± coming from the e + e~ pair. Finally at least 30 GeV energy deposited in one or both of 
the ZDCs was required to select events with forward neutron emission. 

The total number of events collected by the UPC trigger was 8.5 M, of which 6.7 M 
satisfied standard quality assurance criteria. The useable event sample corresponds to an 
integrated luminosity £i nt = 141 ± 12 pib~~ . 

The following cuts were applied to enhance the sample of genuine 7-induced events: 

1. A standard offline vertex cut |wta; z | < 30 cm was required 

2. Only events with exactly two charged particles were analyzed. This cut allows to 
suppress most of non photoproduction contamination in the UPC trigger. 

3. A RICH cut selects which fire 2 or more tubes separated by the nominal ring radius. 

4. Good Track-EMCal matching is also required. 

5. An EMCal energy cut {E\ > 1 GeV || E 2 > 1 GeV) is applied to select candidate e* 1 
above the ERT trigger threshold. 

6. Events with back-to-back e + e~ candidates (detected in opposite arms) were selected. 
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Figure 2: Left: (a) Invariant mass distribution of e + e~ pairs fitted to the combination 
of (shaded) a dielectron continuum [exponential distribution] and (cross hatched) a J/ip 
[Gaussian] signal. The two additional dashed curves indicate the maximum and minimum 
continuum contributions considered in this analysis (see text), (b) J/ip invariant mass 
distribution after subtracting the fitted dielectron continuum signal in (a). 



3 Results 

After the above cuts we find 28 events with e + e~ pairs and none with like-sign pairs for 
m e ± e ± > 2 GeV/c 2 . The measured e + e~ invariant mass distribution for the sample is shown 
in Fig. [5] a) . This distribution is fitted with a continuum (exponential) curve combined with 
a Gaussian function at the J/tp peak, as shown by the solid curve in Fig. [2] a). Simulations 
based on events generated by the starlight Monte Carlo [13l[T7l[T8] and processed through 
GEANT3 [19] have shown that the shape of the measured continuum contribution is well 
described by an exponential function dN/dm e + e - — A ■ e cm " + ^~ . Those simulations allow 
us to fix the exponential slope parameter to c = —1.9 ± 0.1 GeV _1 c 2 . The combined data 
fit is done with three free parameters: the exponential normalisation (A), the J/ip yield 
and the J/tp peak width (the Gaussian peak position has been fixed at the known J/ip 
mass of mj/^ = 3.097 GeV/c 2 [20]). The J/tp and continuum yields and the corresponding- 
statistical errors are calculated from the fit and summarized in Table [T] Fig. [2]b) shows the 
resulting invariant mass distribution obtained by subtracting the fitted exponential curve of 
the dielectron continuum from the total experimental e + e~ pairs distribution. 

Physical cross-sections were obtained after correcting the raw number of signal counts for 
the geometrical acceptance of our detector system, and the efficiency losses introduced by 
the previously described analysis cuts. Acceptance and efficiency corrections were obtained 
with a full Monte Carlo of the experimental apparatus with realistic starlight Monte 
Carlo. Such a model reproduces well the existing d 3 N/dyd<pdpT distribution of coherent 
p production in UPC Au+Au events measured at RHIC by STAR [10]. The coherent 
events were simulated in the PHENIX detector using GEANT3 and passed through the same 
reconstruction program as the real data. 
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m e+e - [GeV/c^] 


Yield 


Cross-section [/ib / (GeV / 'c z )\ 


J/ip peak 


9.9 ±4.1 (st) ± 1.0 (sy) 


76 ± 31 (st) ± 15 (sy) [/ib] 


e+er cont. [2.0,2.8] 


13.7 ±3.7 (st)±1.0 (sy) 


86 ± 23 (st) ± 16 (sy) [starl.: 90] 


e+e- cont. [2.0,2.3] 


7.4 ±2.7 (st)±1.0 (sy) 


129 ±47 (st)±28 (sy) [starl.: 138] 


e+er cont. [2.3,2.8] 


6.2 ±2.5 (st)±1.0 (sy) 


60 ± 24 (st) ± 14 (sy) [starl.: 61] 



Table 1: Measured J/ip and e + e~ continuum photoproduction yields and cross-sections 
at midrapidity in ultra-peripheral Au+Au collisions (accompanied with forward neutron 
emission) at y/s NN = 200 GeV (obtained from the fit of the data to an exponential plus 
Gaussian function) per invariant mass range are shown in the second column. For the 
continuum cross-sections starlight predictions are taken from Ref. [TS] , 



For J/tp at midrapidity the differential cross section is calculated as: 

|y|<0.35 BR Acc ■ E ■ £ tn gg ■ £int Ay 



d(Jj/^+Xn 



dy 



Here Acc is the detector acceptance, e is the track reconstruction efficiency, £trigg is the 
trigger efficiency. The integrated luminosity Ci n t is given in Section 2, Ay is the rapidity 
interval of the measurement. These correction factors and corresponding uncertainties are 
quoted in Ref. [5], and BR = 5.94% is the known J/ip dielectron branching ratio [2"0] , 

For dielectrons at midrapidity (y is the rapidity of the pair) the double differential cross 
section is: 



d 2 a P + 



e~+Xn 



dy dm e 



N e+e - 1 l_ 

|l/|<0.35, Am e+[ _ Acc ■ E ■ Etrigg ■ Cint Ay Am e +, 



where the factors are defined as for the previous equation. 

The measured dielectron cross sections at midrapidity are in very good agreement with 
the STARLIGHT predictions for coherent dielectron photoproduction (rightmost column of 
Table [J) [H]. 

The final J/ip cross section is in good agreement, within the (still large) statistical errors, 
with the theoretical values computed in p~7l [HI [T3l EE HB1 HI] as shown in Fig. [3j For details 
of the comparision sec Ref. [5] 
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Figure 3: Measured cross section of J/ip + Xn production at midrapidity in UPC Au+Au 
collisions at ^s NN = 200 GcV. The error bars (boxes) show the statistical (systematical) 
uncertainties. When available, the theoretical calculations for the coherent and incoherent 
components are shown separately in (a), and summed up in (b). 
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